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In this study the Singlet Reference Interaction Site Model (SRISM) is employed to the study of the electrode 
charge dependence of the capacitance of a planar electric double layer using the primitive model of the 
double layer for a high density electrolyte that mimics an ionic liquid. The ions are represented by charged 
hard spheres and the electrode is a uniformly charged flat surface. The capacitance of this model fluid is 
calculated with the SRISM approach with closures based on the hypernetted chain (HNC) and Kovalenko- 
Hirata (KH) closures and compared with simulations. As long as the magnitude of the electrode charge is not 
too great, the HNC closure shows the most promise. The KH results are reasonably good for a high density 
electrolyte but are poor when applied at low densities. 
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1. Introduction 

In a recent paper [1J, two of us together with Bhuiyan (WBH) applied the Singlet Reference 
Interaction Site Model (SRISM), using both the hypernetted chain (HNC) and Kovalenko-Hirata 
(KH) closures, to the study of density profiles of the double layer formed by a restricted primitive 
model (RPM) at low electrolyte densities and concluded that the HNC closure was preferable. 
However, that study was an introductory one whose purpose was to apply and test our programs 
to the previously studied systems. 

The RPM is a simple, but informative, model of an electrolyte that models the ions as charged 
hard spheres, all of the same diameter, in a dielectric continuum that might represent the solvent 
in the case of ions dissolved in a solvent or represent the environment of the ions or a diminished 
partial charges of the ions in the case of a high density system. It has been employed in the study 
of double layers for aqueous or other solvents. An electrical double layer (EDL) is formed by the 
ions near a charged electrode because the charge of the electrode attracts a compensating layer of 
counterions. This compensating layer may be a single layer or may consist of several layers, some 
of which may consist predominantly of coions. 

In this paper, we study the model electrolyte first introduced by Lamperski et al. [5] and recently 
considered by Outhwaite et al. [3 . This model is based on the RPM and consists of charged hard 
spheres whose diameter is d = 4 A and whose dimensionless temperature is T* — edk^T/ z 2 e 2 = 0.8, 
where fee is the Boltzmann constant, T is the temperature of the system and e is the magnitude 
of the elementary charge. The ions are assumed to be symmetric with charges ze Q and — ze Q and 
e = 78.5 is the background dielectric constant of the environment of the ions. Lamperski et al. 
(LOB) examined the capacitance of this model system by means of a simulation and found that 
the capacitance had a minimum at zero electrode charge at small electrolyte densities. Such a 
minimum is typical of the Gouy-Chapman-Stern (CCS) theory [IHB], that neglects ion diameters, 
except for the distance of closest approach of the ions, and correlations between the ions. This 



© S. Woelki, D. Henderson, 2011 



43801-1 



S. Woelki, D. Henderson 



parabolic-like minimum is seen in simple low density electrolytes. As the density is increased, this 
minimum becomes shallower but is always present in the GCS theory. By contrast, LOB found, 
for their simulations that used their model system, that at higher densities the minimum in the 
capacitance fills and becomes a maximum. As pointed out by Kornyshev [7], such a maximum is 
typical of many ionic liquids that exhibit a bell shaped or dromedary camel- like capacitance curve. 
Kornyshev's theoretical observations have been supported by recent experimental work by Lockett 
et al. [5] and Islam et al. [5J[TU]. Indeed, the LOB model was not meant to represent any specific 
system but docs mimic an ionic liquid. LOB also found that the modified Poisson-Boltzmann 
(MPB) theory [TT] exhibited this same behavior. Since the MPB theory fails to converge at high 
electrode charge, they did not examine the capacitance at high electrode charge. Subsequently, 
Lamperski and Henderson (LH) [T^] extended this study to a very large electrode charge and found 
in their simulations that the capacitance at low densities has a minimum at zero electrode charge 
but the capacitance has a maximum at moderate electrode charge. This maximum is followed by 
a steadily decreasing capacitance. The capacitance at low densities has the shape of the letter 'm', 
or a bird with drooping wings, or a bactrian camel. At large electrode charge, the capacitance in 
the LH study is density dependent but is not strongly so. Thus, at all electrolyte densities, the 
capacitance decreases steadily at a high electrode charge. This is due to a thickening of the double 
layer because there is a limit to the amount of counterion charge that can be accommodated near 
the electrode. By contrast, the GCS theory neglects the ion size and correlations and there is no 
limit to the amount of charge that can approach the electrode in this theory. 

The purpose of this paper is to apply the SRISM theory of our previous paper to the LOB 
system. The present study is an interim stage in our work since SRISM is not restricted to spherical 
ions. Indeed, our long-range intention is to study double layers that contain nonspherical ions but 
our plan is to proceed in a step by step manner. 

2. Theory 

We reviewed the SRISM theory in detail in our previous paper [T] and there is little need to 
repeat this discussion. Two closures were examined, the hypernetted chain (HNC) and Kovalenko- 
Hirata (KH) closures. The KH approximation is a hybrid of the HNC approximation and the 
mean spherical approximation (MSA). These closures employ the Ornstein-Zernike (OZ) relations, 
given by equation (2) of WBH, that connects the total and direct correlation functions. The total 
correlation functions are the radial distribution functions, less their asymptotic values (unity). 
Thus, 

hij(r) = gij(r) - 1, (1) 

where the hij(r) are the total correlation functions and the gij(r) are the radial distribution 
functions that give the probability of finding ions of species i at a distance r from a central 
ion of species j. The direct correlation functions, Cjj(r), are defined by the OZ relations. As the 
name implies, the total correlation functions give a complete information about all the correlations 
between a pair of ions whereas the direct correlation functions give the correlations of such a pair 
directly, but in the presence of the other ions. Since the OZ relations are a definition of the direct 
correlation functions, they do not yield a method of calculating any of these correlation functions. 
The OZ relations must be supplemented by an approximate closure. For example, for the spherical 
system considered here, the HNC closure is 

hij(r) = dj{r) + lnpy(r) + fiuij(r), (2) 

where f3 — with fee being the Boltzmann constant and Uij{r) is the pair potential between 

a pair of ions of species i and j and charge Zie and Zje a . The pair potential is given by 

««(r) = ^ for r^d 

er 

u ij( r ) — 00 f° r r < d. (3) 
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The KH closure is a hybrid of the HNC closure and the MSA closure, which is a linearized 
version of the HNC closure, and is given by 

Cij (r) = — ftiiij (r) for r ^ d (4a) 

hij(r)=-l for r < d. (4b) 

The second part of the above equation is merely a statement of the fact that the ions are impene- 
trable. The first part is the approximation. Specifically, the KH closure is 

Cy(r) = - @Uij(r) for hy (r) - cy (r) - /3ity (r) ^ (5a) 

Cij(r) = hi j(r) -In gij(r) - /3uij(r) for fty (r) - cy (r) - (r) < (5b) 
hij(r) = — 1 for r < d. (5c) 

Again, the third part of the above equation is only a statement of the fact that the ions are 
impenetrable. 

Our procedure is to calculate the pair correlation functions of the bulk electrolyte by solving 
the bulk OZ equations with the chosen closure and then using these results as an input to solve 
the surface OZ equations for the singlet surface correlation functions. The ion-surface interaction 
potentials are 

u±(x) = — ^ 7TZ±e ° a x f or x ^ 

e 

U±{f£)=OQ for x < d/2, (6) 

where x is the normal distance of the center of an ion from the electrode and a is the electrode 
charge density. In some applications there is an advantage in using a separate closure for the 
computation of the bulk electrolyte correlation functions. This occurs when the MSA is used for 
a bulk electrolyte with spherical ions because analytical results can be used. We will not do this 
here and instead use the same closure for the bulk and interfacial parts of our calculation. 

Once the interfacial correlation functions have been determined, the interfacial charge density 
and electric potential drop across the double layer can be determined by integration. For example, 



a = —4ire p J [h + (x) — h- (x)} dx (7) 

x=0 



and 

4ire a p 



p / [h+ (x) - h- (x)] xdx. (8) 



x=0 

We specify the charge density a as the input variable so equation ^ is merely a check. Once <f> as 
a function of a has been determined, the integral capacitance, 

C,= a -, (9) 

and differential capacitance, 

cr d = |, (io) 

can be obtained. We use the dimensionless electrode charge density, a* = ad 2 /e and electrode 
potential, 4>* = f3e <fi. Following LH, we report values for the dimensionless differential capacitance, 
C d d. 
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3. Results 

In figure [T] we show our results for differential capacitances for the three densities, p* = 
Nd 3 /V = 0.04,0.14,0.24, considered by LH. Additionally, we show our HNC results for a higher 
density p* = 0.34 in figure |] By MSA/MSA1 we mean that the MSA is used for the bulk corre- 
lation functions and the singlet interfacial correlation functions, respectively. The SRISM theories 
fail for large a* so we restrict our comparison to a* between —0.3 and 0.3. The MSA results are 
quite reasonable at very low electrode charge. This has already been reported by some of us [T3] • 
However, because the MSA is a linear response theory, the MSA capacitances are independent of 
a*. The KH results are very poor for the p* = 0.04 system. They are rather like the MSA re- 




-0,3 -0,2 -0,1 0,1 0,2 0,3 

a* 

Figure 1. (Color online) The reduced differential capacitance, = Cdd, of the electrical dou- 
ble layer as a function of the reduced surface charge density, a* for d = 4 A at the reduced 
temperature T* = 0.8 and the reduced electrolyte densities (a) p* = 0.24, (b) p* = 0.14 and (c) 
p* = 0.04. The circles are the simulation results from reference |12| and the black (middle), red 
(upper) and blue (lower) lines show the corresponding MSA/MSA1, HNC/HNC1 and KH/KH1 
results, respectively. 

suits because they change little with a* . However, the KH results improve with increasing density. 
Indeed, the KH closure correctly predicts a maximum in the capacitance for p* — 0.24. Further 
study is needed to determine the general utility of the KH closure for high density (ionic liquid) 
systems. The HNC results are fairly reasonable for p* = 0.04. For the highest densities, the HNC 
predicts only a very weak maximum at p* — 0.24. However, it is at least a maximum. As is seen 
in figure [2j there is a clear maximum in the HNC capacitance at the highest density. Simulation 
results for this high density are not available and probably unnecessary as the purpose of figure [2] 
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Figure 2. The reduced differential capacitance, = Cdd, of the electrical double layer as a 
function of the reduced surface charge density, a* for d — 4 A at the reduced temperature 
T* = 0.8 and the reduced electrolyte density p* = 0.34. The black (lower) and red (upper) lines 
show the MSA/MSA1 and HNC/HNC1 results, respectively. 

is merely to demonstrate that the HNC theory does exhibit a maximum in the capacitance at a 
sufficiently high density. The relevant point is that the HNC theory predicts a maximum at high 
densities but it is delayed. 

4. Summary 

We have applied the SRISM approach to the study of the capacitance of the double layer formed 
by charged hard spheres, all being of the same diameter and magnitude of charge, near a flat 
uniformly charged unpolarizable surface for a range of densities and electrode charge. Comparison 
with simulations shows that the KH is rather poor at low electrolyte densities but improves as 
the density is increased. The HNC is fairly good overall but the appearance of a maximum in 
the capacitance is delayed and does not appear clearly until rather high densities. The MSA is 
fairly good for a small electrode charge but is incapable of accounting for the electrode charge 
dependence of the capacitance due to of the linearization inherent in the approximation. None of 
the theories considered here is capable of predicting the decreasing values of the capacitance at a 
large electrode charge because they neglect electrode correlations and permit too many ions near 
the electrode and predict a double layer that is too thin. 

On the whole, the MPB results are somewhat better than those reported here. However, since 
the MPB theory, at least in its present incarnation, fails to converge when the electrode charge is 
large in magnitude, it is not known whether the MPB theory predicts the continuing decrease in 
the capacitance at large electrode charge. Also, at present, the MPB theory provides results only 
for spherical ions. The results presented here are not unpromising. The HNC approximation, at 
least, seems useful as long as the electrode charge is not too large. The KH closure may be useful 
at high electrolyte densities although a further study is needed to determine whether this is indeed 
the case. Finally, the SRISM approach can be applied to nonspherical ions. 

The recent work by Loth et al. |14| presents an alternative theory and some simulations for the 
capacitance of a high density fluid. 
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3acTOcyBaHH5i MeTOfly SRISM p,o bi/ibhghhji e/ieKTpoMicTKOCTi 
noflBiwHoro uiapy b h co ko rycT m h h oT npuMiTi/iBHOi' 
MOfle/ii e/ieKTpo/iiTy 

P,. TeHflepcoH =1 

I I HCTHTyT 4>i3HHHoT i TeopeTi/i h hoi ximm, y h iBepcuTeT m. PereHc6ypra, HiMeHHHHa 
cfaKy/ibTeT xiMiT i 6ioxiMiT, yHiBepcwTeT Bpirewia 9\nra, ripoBO, CLUA 

B LibOMy flocnifl>KeHHi MOfle/ib cm Hr/ieTHoro BUflijieHoro B3aeMOfliKDHoro LieHTpa (SRISM) 3acTocoBye~rbca flo 
BMBHeHHiR 3ajie>KHOCTi e/ieKTpoMicTKOCTi n/iocKoro e/ieKTpuHHoro noflBiMHoro inapy Bifl 3apafly e/ieKTpofla, 
BUKopucTOByioHi/i n pi/iMiTi/iBHy MOfle/ib noflBiwHoro ujapy fljia Bi/icoKorycTHH Horo ejieKTpo/iiTy, iu,o MOflejiioe 
ioHHy pifli/iHy. Iohm npeflCTaBJiaioTbca sk 3apafl>KeHi TBepfli ccpepw, a e/ieKTpofl e oflHopiflHO 3ap:qfl>KeHoio 

n/lOCKOKD nOBepXHeK). MicTKicTb LlbOTO MOflejlbHOrO njlUHy p03paXOByeTbC51 MeTOflOM SRISM i3 3aMHKaHHiRMI/1, 

u^o 6a3yKDTbca Ha rinep/ia Hi4KD>KKOBOMy (HNC) niflxofli i 3a m n Ka h H3X KoBa/ieHKa-ripaTH (KH) i nopiBHtoerbOR 
3 cwMy/iflLiiflMH. 51klho Be/iHHWHa 3apafly ejieKTpofla He e flyx<e Be/i n kohd, HNC 3a MHKa h H3 BMnnaflae HaR6i/ibiu 
nepcneKTi/iBHi/iM. KH pe3y/ibTaTH e npwtiHflTHi/iMi/i p,nn e/ieKTpo/iiTy npi/i bwcokmx rycTi/iHax, ajie ripiuwMH npi/i 
Hi/i3bKwx rycTHHax. 

K/iiOMOBi c/ioBa: cnHrneTHi 'mTerpanbHi piBHXHHn, ioHHa piflviHa, oBMemeHa nptAMiTtABHa MOftenb, 

ejieKTpOMicTKICTb 



43801-6 



